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Time Evolution of the Large-Scale Tail of Non-Helical Primordial Magnetic Fields 
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We present a derivation of the time evolution equations for the energy content of non-helical 
magnetic fields and the accompanying turbulent flows from first principles of magnetohydrodynamics 
in the general framework of homogeneous and isotropic turbulence. This is then applied to the Early 
Universe, i.e. the evolution of primordial magnetic fields. Numerically integrating the equations 
we find that most of the energy is concentrated at an integral wave number scale kj where the 
turbulence turn over time equals the Hubble time. At larger length scales L, or smaller wavenumbers 
q = 2n/L <C ki, independent of the assumed turbulent flow power spectrum, mode-mode coupling 
tends to develop a small q magnetic field tail from zero with a violet noise slope proportional to the 
fourth inverse power of L and therefore a scaling for the magnetic field of B ~ L~ 

PACS numbers: 95.30.Qd, 98.62.En, 98.80.Cq 
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I. INTRODUCTION 

The question of the origin and time evolution of pri- 
mordial magnetic fields in the Early Universe is an in- 
teresting and so far at best partially resolved problem in 
cosmology. It is possible that strong magnetic fields have 
been created on small scales in the Early Universe by a 
cosmological process, for example at a phase transition 
[l[ or, less likely, during inflation Q (for an overview of 
ossible magnetogenesis scenarios see, for example, Q or 
). A central question then is whether these magnetic 
fields could evolve with time and transport some of their 
energy content to large scales to account for the recently 
claimed detection of Intergalactic Magnetic Fields 0, Q . 

Many attempts to study the evolution of primordial 
magnetic fields have been performed in the past 043- 
Numerical simulations are problematic as they lack the 
resolution required to give reliable predictions. A further 
complication in the study of the evolution of magnetic 
fields in the Early Universe is the enormous cosmic ex- 
pansion between a putative magnetogenesis scenario and 
the present, such that the smallest errors in extrapolation 
lead to large changes in the final prediction. 

In this paper we take a semi-analytic approach which 
derives the main time evolution equations from first 
principles of magnetohydrodynamics (MHD), employing 
some fairly generic assumptions. Our analysis follows a 
similar procedure as has been already applied to, for ex- 
ample, the solar wind [l7|, [U| or the galactic dynamo 
problem pi]. 

Two of us [2(| have recently, for the first time, applied 
such an approach to the evolution of primordial mag- 
netic fields. Considering only the most important large q 
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velocity-magnetic field mode-mode coupling source term 
to generate magnetic fields on small q taken from (l9j . 
and building on the result of , it has been established 
that cosmic expansion seems slow enough to allow for 
the generation of large scale magnetic fields with a white 
noise, i.e. B ~ L~ 3 / 2 , or even shallower spectrum, also 
in the absence of initial magnetic and velocity fields on 
such large scales [2l[. However, the analysis is far from 
complete as the role of back-reaction of the turbulent 
medium onto the magnetic fields has not been accounted 
for. Here we present the full analysis including all source 
and sink terms for non- helical fields. 

This work is structured as follows: In Section [H] we 
give the evolution equations for the magnetic and tur- 
bulence fields in the general framework of homogeneous 
and isotropic non-compressible MHD. An outline of the 
lengthy derivation of these equations is delegated to Ap- 
pendix [AJ whereas a different form more suitable for nu- 
merical integrations is presented in Appendix [Bj In Sec- 
tion UTT] this is then extended to the situation in the Early 
Universe with the results and conclusions from these con- 
siderations finally presented in Sections HVl andlVl respec- 
tively. 



II. TIME EVOLUTION OF THE MAGNETIC 
AND KINETIC ENERGY CONTENT IN 
HOMOGENEOUS ISOTROPIC 
MAGNETOHYDRODYNAMICS 

The equations describing the time evolution of the two 
main observables of magnetohydrodynamics, the velocity 
field of the turbulent medium v and the magnetic field 
B, are, for an incompressible fluid (i.e. Vv = 0), given 
by 



d t B 



1 

47TCT 



AB + V x (v x B) 



(1) 



2 



and 

a t v = -(vV)v+ (V * B)xB + f,, (2) 

respectively. Here a is the conductivity (which is as- 
sumed to be very large in the following) and p the mass 
density of the fluid, while f v is some viscous density force. 
Note that for most part of the evolution of the early Uni- 
verse incompressibility is an excellent assumption due to 
the large speed of sound in the relativistic plasma. 

Our main interest lies in the average buildup of mag- 
netic and kinetic energy density, respectively, taken over 
an ensemble of cosmic realizations (denoted by chevrons 
()), i.e. 

(d t M q ), {d t U q ) (3) 

with M q being the magnetic spectral energy defined 
through 

£B = 8^/ d3 " B2(x)= / ^|B(k)| 2 ^p/dfcM fe 

(4) 



and U q being the kinetic spectral energy given by 

€K = WJ dW ( x ) = f / d 3 fc|v(k)| 2 ^p J dkU k . 

(5) 

For both expressions V denotes the volume and we have 
assumed cosmic homogeneity and isotropy which implies 
that Mfc and U k are functions only of the magnitude k of 
the wavevector k. Furthermore, Parseval's Theorem has 
been used in both cases in order to obtain a k integral 
where a hat denotes the Fourier Transform normalized 
by V? (cf. Appendix lAl . With these assumptions the 
spectral energy densities are therefore 

M g = |^|B(q)| 2 (6) 

and 

U q = 2izq 2 \v(cj)\ % . (7) 

By performing the calculations which are presented in 
Appendix |Al we obtain a very general result for the ho- 
mogeneous and isotropic case which for ([3]) is given by 



(d t M q ) 



dki At / d9 



i) 



sin 3 9 (M q ) (U kl ) + \^(q 2 + k 2 -qk cos 9) sin 3 9 (M k ) (U kl ) 

Ivy rZ-y 



-q 2 (3 - cos 2 9) sinO (M k ) (M q ) 



(8) 



and 



dki At / d6 



- § (qk sin 2 9 + 2k\ cos 9) sin 9 (M k ) (M kl ) + - ^ (3fc - q cos 9) sin 3 9 (U k ) (U kl } 

4 fvy *-± rZ-y 



1 qk 4 



1 k 4 
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4 V ^ 3<1 + k 9) Sin3 6 {Uq) {Ukl ^ ~ 2 kf ^ + ^ ~ ^ C ° S ^ Sin3 ° {Mkl ) '' Uq)+ 2 k2 Sin3 6 (Mfe) {Uq) 



(9) 



where q, k and k\ are the magnitudes of the wavevectors 
q, k and ki = q k, respectively, and 9 is the angle be- 
tween these two vectors, i.e. qk = qk cos 9. Eqs. ([8]) and 
© are a set of well-defined equations since they ensure 
conservation of energy, momentum and mass density to 
the lowest non-trivial order in At [13] ■ 



III. APPLICATIONS TO THE EARLY 
UNIVERSE 

In order to study the time evolution of magnetic fields 
in the Early Universe we include Expansion. For the fol- 
lowing considerations it is convenient to introduce the 



scale factor a as the time evolution parameter which in 
this work is set to unity at the initial magnetogenesis era. 
The focus here will be on the long radiation dominated 
period in the Early Universe where most of the non-trivial 
evolution occurs. The derivation of the MHD equations 
in an expanding Universe has, for example, been pre- 
sented in [9(. These equations contain redshifting terms 
proportional to powers of the scale factor, properly ac- 
counting for, e.g., the physical decrease of magnetic fields, 
i.e. B ~ 1/a 2 due to flux-freezing. These terms are not 
taken into account in Eqs. ([8]) and ([9]) and have to be 
included now. However, it is known that with a proper 
scaling of variables (cf. Appendix B in [3), 



3 



dt r = dt ( 



- T3 „2 



B r = B 
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pa k c = ka , 
. ( 10 ) 

the MHD equations in an expanding radiation dominated 
Universe are forminvariant to the corresponding equa- 
tions in a static background. Our master equation is 
therefore also valid in the Early Universe when these " co- 
moving" (denoted above by a superscript "c") variables 
are used. It is convenient to change the time derivative 
to a scale factor derivative 



d_ 

~dt c 
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2toa d In a 



h d 

a d In a 



(11) 



where to and Hq are cosmic time and Hubble parame- 
ter at the initial magnetogenesis period, respectively. It 
is noted here that due to our definition of the Fourier 
Transform given in Appendix [Al the comoving Fourier 
Transforms for B c and v c are given by 



B' 



Ba 1 / 2 



-3/2 



(12) 



This implies that M c 



M q a 



,,, ,^ 1 and Ug = Uqa^ 1 , such 

that e c B and e c K are constant during the expansion of the 
Universe, when dynamical evolution is excluded. In what 
follows, the superscript "c" is dropped and all variables 
are comoving, unless specifically noted otherwise. 

It is now possible to rewrite (J5J and © in comoving 
coordinates as 



d In a 

dU q 
d In a 



Ho 
a 



f ik { At I 



d6 



(13) 



(14) 



where the square brackets [...] in (|T3"|) and (|T4")) denote 
remaining terms in (|8]) and ([9]), respectively. 

Finally, a choice for At has to be made. Following 
methods applied in molecular chaos we choose At to 
be the smallest of all "eddy" turnover times, i.e. At ~ 
min(L/u L , L/va) where v L and v\ are the effective fluid 
and Alfven velocities on the scale L. In particular we 
take 



At 



2tt 



2tt 



Ho k(2k{U k ))i' k(3/2k(M k ))\ 



, (15) 



energy including back-reaction by numerically integrat- 
ing Eqs. (fT3"j) and (fTU). We assume two different initial 
conditions. In all cases it is assumed that magnetic fields 
and turbulent velocities are created on some small scale, 
called the initial integral scale k J = ko, i.e. for q = ko 
both (M q ) and (U q ) have a sharp peak. The two cases 
distinguish themselves in the initial conditions for the 
small-g tail of the turbulent flows, a tail qll q ~ q 5 and no 
U q tail at all. Results of the evolution of M q and U q for 
these initial conditions are shown in Fig. [TJ 

It is striking to note that (a) the peak and the inte- 
gral scale follows the analytical prediction very well, (b) 
the approximate equipartition between magnetic and tur- 
bulent energy on all scales is achieved and (c) the final 
result is independent of the initial conditions, predicting 
a qM q ~ q 5 large scale tail for the magnetic fields. Here 
the last result holds even for the case where a qU q ~ q 3 , 
i.e. an initial white noise turbulent spectrum is assumed. 
Thus the conclusion concerning the small q magnetic tail 
seems very robust. Below analytical arguments are given 
to support all these results. 



Evolution of the integral scale: 

The simulations follow very well the proposed evolu- 
tion of the peak (i.e. " integral" ) scale kj (a) with scale 
factor a. It has been shown [3, [2(| that if the spec- 
tral energy has a slope proportional to g Q_1 {a > 1) for 
q < k], i.e. 



(E q ) 



(16) 



where E can be either M or U, and Eq (i.e. Mo or Uq) 
is the normalization constant given by the corresponding 
initial value at ko, we have 



ki « koa a + 2 



(17) 



for the time dependent integral scale, ko being deter- 
mined by [20l | 



ko 



2itH 



i'o 



27T 2 if 2 

U 



(18) 



and 



such that we identify (2k (U k )) 2 and (3/2fc(M fc )) 2 with 
the effective fluid and Alfven velocities on scale L — 
2ir/k. Note that the first timescale in the brackets of 
Eq. fp~5|) is supposed to ensure causality, in practice, how- 
ever, our results are independent of this condition. 



E a 2 °+ 2 



(19) 



for the value of the corresponding spectral energy at kj. 
Here, for initial equipartition with vq ~ 1, it has been 
assumed that E = M = Uo = (2fc ) _1 . 



IV. RESULTS 



Large-scale magnetic tail: 



With this at hand it is now possible to analyze the time 
development of both the spectral magnetic and kinetic 



Large scale (i.e. small q) tails are best discussed in 
terms of Eqs. (|B1[) and (|B2[) . For large scales, i.e. q <^L 
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FIG. 1. Time evolution of magnetic (black) and kinetic (grey) 
spectral energies according to ()13|) and (|14[) . respectively. 
Dashed lines denote the initial conditions while the solid lines 
represent the situation for a = 10 s . Upper panel: Starting 
at some time where both spectral energies were concentrated 
on the same scale they evolve close to equipartition, building 
up a q 4 slope ab initio. Lower panel: Starting at some time 
when the turbulence has already built up an q 4 slope after 
some time the magnetic spectral energy reaches equipartition 
with the same slope fairly well. It should be noted that for 
both cases the actual value for ki from the simulation is in 
good agreement with the one predicted by (JT7J), denoted by 
the vertical lines labeled 'initial ki' and 'final ki\ respectively. 



kj ~ k, we can assume (U q ) <C (Uk) and (M q ) <C (Mk) 
and therefore neglect the terms containing (M q ) and (U q ) 
in (|BT|) and (|B2|) . leaving the (M k ) (U kl ) term in the 
former and both the (Mk) (M kl ) and (Uk) (U kl ) terms in 
the latter one. 

If we now use the power law assumption motivated 
above, i.e. (M q ) ~ q a ~ x and (U q ) ~ q@~ l , it is possible 
to evaluate the k\ integrals which, performing a Taylor 
series in q/k, give 



' 8M q ' 
9 In a , 



3 (M k ) (Uk) 



O 1 



dU q 
d In a 



dfc F 



o 



(Uk) (Uk) - ^ (M k ) (M k ) 



(21) 



One of the main results from these considerations is the 
fact that for q <C fc/ both the magnetic and the kinetic 
spectral energies have a steep spectrum which is propor- 
tional to q 4 or L~ 4 , i.e. a — (3 = 5, where L is the scale 
of the corresponding field. This would correspond to a 
behavior of the form B ~ L _ 2 and v ~ for the scale 
dependence of the magnetic and the turbulence fields on 
large scales, respectively. For the former this would be in 
good agreement with [23[ and [24{ which predict the same 
scaling. For the turbulent field it is an important result 
which is also confirmed by simulation as it can be seen in 
Fig. [T] Fairly independent of the initial condition, even 
starting off with the total energy being concentrated on 
only one scale, a q 4 slope for the spectral energies forms. 



Equipartition: 

Furthermore, the system tends to achieve equipartition 
quite accurately. A simple argument for this at large 
scales can be extracted from (|2TIt and ([2"Tj) : Assuming 
that the magnetic and kinetic spectral energies have sim- 
ilar values at k, i.e. (Mk) ~ (Uk), and a — 5 as discussed 
before, we obtain 



dU q 
d In a 



dM q 
d In a 



dk 



{Uk) {Uk) 



which means that the time evolution of both is the same 
to the fourth order in q. 

This is a crucial result since it means that equipar- 
tition is indeed possible or at least that at some point, 
even for the case that M q and U q do not have the same 
magnitude, the time development of both at large scales 
will qualitatively be the same. 



Estimate of present magnetic field strength: 

To estimate the magnetic field strength B(L) at some 
scale L = 2-k jq we use 



M„ 



B\L) 



(20) 



8irqp 

and, with ([T7|). (ITS)) and ([15]). obtain 

B(L) = (8TrqpM q )i = B (2qM q )i w B (H L)~ 



(23) 



(24) 



5 



where Bo = (47rp) 2 is the effective magnetic field for vq « 
1, i.e. the speed of light, for the initial equipartition of 
radiation and magnetic energies which would correspond 
to So ~ 3 x 1CP 6 G. If we now acknowledge the integral 
scale to be the coherence scale of the magnetic field, then 
for the magnetic field strength at Lj = 2n/ki we get, 
using (fTT ]) and (fT9|) . 



A(x) by 



B(Lj) = B (2k I M I ) 2 = Boa^+i 



(25) 



and therefore for a — 10 8 and a = 5 it is Lj w 200 pc 
and hence B(200pc) < 5 x 10~ 12 G for the QCD phase 
transition. It is hard to imagine how causally generated 
magnetic fields in the early Universe could yield a larger 
remnant field than quoted above, unless they are gener- 
ated with substantial helicity. 



V. CONCLUSIONS 

Starting from first principles and reasonable as- 
sumptions about homogeneous and isotropic magneto- 
hydrodynamic turbulence, evolution equations for the 
spectral magnetic- and kinetic- energy densities have 
been derived. After adapting these equations to the 
use in an expanding Universe they were numerically in- 
tegrated for a number of assumed initial conditions, as 
possibly resulting during an early magnetogenesis period. 
It has been found that, seemingly independent of initial 
conditions, a B(L) ~ L -5 / 2 tail on scales L ^> 2n/ki al- 
ways develops. This magnetic field spectrum may there- 
fore be regarded as the natural one for non-helical cosmic 
magnetic fields (cf. also to (24J). At the same time most 
of the energy is concentrated on the integral scale which 
allows a rather generic prediction of coherence scale and 
strength of causally created primordial magnetic fields. 
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Appendix A: Derivation of the Master Equations ([8]) 
and © 

The first step in order to derive © and © is to obtain 
an expression for ^ and ([7]). To do so we first take the 
Fourier Transforms of the initial equations (JXJ) and ©, 
where we define the Fourier Transform of a vector field 



A(x) 



V* 
(2n)i 



d 3 fcA(k)e 



kx 



(Al) 



Using the convolution theorem, for a wavenumber q they 
are given by 



Wq) 

and 
Stv(q) 



1 

Ana 
+ (2tt)§ 



<Z 2 B(q) 



d 3 fc|q x v(q-k) x B(k) j 



(A2) 



+ 



iV? 1 
(2tt)* 4^ 



d 3 /c{ k-v(q-k) v(k)j 
d 3 fc{ 



kxB(k) xB(q-k)}. 

(A3) 



We solve these ordinary differential equations by using 
the Midpoint Method which gives us B(q, t) and v(q, t) 
at the time t = to+At for some initial conditions Bo(q) = 
B(q, to) and vo(q) = v(q, to). In particular, if the time 
derivative of some function A is dtA = f(t,A(t)), the 
Midpoint Method is given by 

A(t+At)=A(t)+Atf (t+At/2, A(t)+At/2 f t, A(t) 

We use this method as it is accurate to second order 
in the time step At, essential for energy conservation of 
Eqs. (0) and ©. 

In order to calculate §5§ we use ^ and ([7]) and there- 
fore obtain 



(d t M q ) 



2p 



|B(q)| 2 



(£_ 

2p 



9 t B(q) ■B(q)*)+h.c. 



q 2 B(q, t) ■ B(q, t)* - B(q, to) ■ B(q, i )' 



2p 



At 



and 

(dtU q ) 

= (d t [2V|v(q)| 2 ]) = (2V[9 t v(q)] • v(q)' 
2 v(q, t) ■ v(q, t)* - v(q, t ) ■ v(q, t )* 



2nq z 



At 



(A5) 



h.c. 



(A6) 



which are accurate to order At if the Midpoint Method 
is used. 

When the terms B(q, t) ■ B(q, t)* and v(q, t) ■ v(q, t)* 
with t = to + At in Eqs. (|A5I) and (IA6|) are evaluated 
with the help of Eqs. (IA2[) and (IA3|) a large number of 
expressions of the sort, for example, as 



6 



(At) 2 ( q 2 



d 3 k / d 3 k' 



q • B(k, t) q ■ B(k', i) v(q - k, f) • v(q - k', t) 



result. Making the reasonable assumption of approxi- 
mately uncorrelated statistical chaos in homogeneous and 
isotropic non-helical MHD turbulence, one may write (cf. 



B a (^t)B b {V ,t')* ) « Ci5w5 w M k S ab - 



(« (k,f)^(k',f)*) « Ci&kk'SwUk Sab ~ 



B B (k,t)f5 i (k',0*) wO 



7) 



(A9) 



where the <Ts denote delta functions and a and 6 are 
constants to be determined. Using Wick's Theorem for 
Gaussian fields, i.e. 

(A Q j4;,j4 c j4d) = (A a Ah) (AcAd) + other permutations 

(A10) 



one may use the relations (|A7j) - (|A9j) to work of one 
integral in expressions of the type of Eq. (|A7j) . yielding 
after computing many terms of this sort the final result 
Eqs. © and ©. Finally it is noted that Ci = p/k 2 
and C2 = p/(47r/c 2 ) as determined form the definitions in 
Eqs. ©and 10}. 



Appendix B: Alternative form of the Master 
Equations 



We give here an alternative formulation of the master 
equations (0) and (0) in terms of ki = q — k which is 
more suitable for numerical integration: 



(d t M q ) 



dki At 



fci, 

kl, m 



d*i 



k 3 



16k 3 k 3 



16k 3 kx 

3 1.3 



/ k\ qk 3 fcf 
^ 16<?fc 3 ~ Sfc 3 " 



9 3 fc 3q 3 



qk 3 
~ Ikf 

q 3 h 



kki 
qk 3 



16kkf 

q 3 h 



16k 3 
5qk 



8qk 16k 3 



i kk-i \ 
8k~ + Wq) < M «> {Mk) 



qk 

Ak[ ' 
qk 



0) (M q )(U kl ) 



8kki 16k 3 16kl 16fci 



qk\ 
16k 



qk 3 \ 
16k 3 J 



(M k )(U kl ) 



(Bl) 



and 

(dtU q ) = 



32k 3 k 3 
fc 7 
32q 3 k 3 

k 7 

16q 3 k 3 
8« 3 fc 







dfcjAt 


/ dfc 




fel.min 


7g 5 


3g 5 


32 kk-^ 


32/c 3 fci 4 


7k 5 


3fc 5 


+ 32qkf 


+ 32 9 3 fci 


k 5 


fc 5 



g 3 /c 3<? 3 



16fcfc 3 

llq 3 k 5q 3 
~32kJ + 
llgfc 3 



8k 3 



3q 3 fc 



gA; 3 
8fcfci ~ 16k 3 

5qk 3 



3qk 
9qk 



16kh 
5k 3 



32k 3 
3k 3 h 



16qk 3 16q 3 ki 



32k\ 32/ci 

5q 3 k 9qk 

' 32kf ~ lQqki ~ 32q 3 + 32k 3 ~ 32h + 

qk 3 3k 3 k 3 ki qk q 3 k kki 
_ + __ _ _ 



5qki 
' "lfjfc 

3qk\ 

~32k ~ 32fc 3 J 
t^kkj^ kk-^ 



(M k )(M kl ) 
qkT] 



j4_ 

Iqk 



kk\ 



qki 

8k 



16kf 
kki 



8qk x 
k 3 h 



16q 3 

(M k ) (U q ) 



32q 

kk\ 
~16q^ 



32q 3 

(M kl 



(U k )(U kl ) 
(U q ) (U kl ) 
(U q )} 



(B2) 



r 



where k\ is the magnitude of ki. The limits of the k\ 
integration, fci >m i„ and ki max , are given by 



fa,min = (q 2 + k 2 - 2qk) 2 , ki >max = (q 2 + k 2 + 2qk) 



(B3) 
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